o 



> 



o 



X 



Hamiltonian analysis of the double null 2+2 
decomposition of General Relativity expressed in 
terms of self-dual bivectors 



o 
o 

CN ■ School of Mathematics, University of Southampton, Southampton, S017 IBJ, UK 

U 

Oh: 

< 



R A d'Inverno, P Lambert and J A Vickers 



E-maih J . A . VickersOmaths .soton.ac.uk 



Abstract. In this paper we obtain a 2+2 double null Hamiltonian description of 

^S| ' General Relativity using only the (complex) SO (3) connection and the components 

of the complex densitised self-dual bivectors Sa- We carry out the general canonical 

analysis of this system and obtain the first class constraint algebra entirely in terms 

tH- ' of the self-dual variables. The first class algebra forms a Lie algebra and all the first 

OO . class constraints have a simple geometrical interpretation. 

O 

o 

\o 
o 



Submitted to: Institute of Physics Publishing 

O^' Class. Quantum Grav. 

U 

bX)' PACS numbers: 0420Fy 



2+2 Ashtekar variables II 2 

1. Introduction 

In two earlier papers [I] (which we will refer to as papers I and II respectively) we 
gave a 2+2 Hamiltonian description of General Relativity using Ashtekar type variables 
based on a self-dual action [Sj, |1I- The variables used in paper II were based on previous 
work by the authors P and also on the approach taken by Goldberg et al [Sj and used 
a mixture of frame variables and the components of the self-dual bivector. This had 
the advantage that the variables had a clear geometric interpretation but suffered from 
the disadvantage that not all the variables were manifestly self-dual. In this paper we 
will describe an alternative approach in which the variables are the components of the 
densitised self-dual bivectors S^ together with the components of the complex SO (3) 
connection 1-forms F"^. As in paper II we use the original Ashtekar formulation in which 
the manifold is real but we allow complex solutions of the field equations. However once 
the canonical analysis has been carried out reality conditions are imposed to limit the 
solutions to real solutions of Einstein's equations. 

In the first section we use the self-dual action introduced in paper I to obtain a 
Lagrangian description expressed in terms of these variables. However if the the bivector 
variables are obtained from a null frame the components are not independent and one 
needs to impose appropriate constraints in order to obtain General Relativity. These are 
found and introduced into the Lagrangian using Lagrange multipliers. This immediately 
gives rise to the primary Hamiltonian and in the next section we look at the Hamiltonian 
description. 

In section 4 we show that both the structure equations for the connection and 
all the Einstein equations are given by the constraints and equations of motion. We 
then go on in section 5 to calculate the first class constraints and obtain the first class 
algebra. This is shown to be a Lie algebra generated by the momentum constraint and 
the modified Gauss constraint. We end by giving a geometrical interpretation to the 
constraints and discuss the next steps in the canonical quantisation process. 

2. 2-|-2 Connection Variables 

As in paper I we we start with the first order self-dual action given by Jacobson and 
Smolin jH] 

'R'^AS^gAB, (2.1) 

where S"^ are the complex self-dual 2-forms, R^ is the curvature 2-form of the complex 
5*0(3) connection F"^ and g^B is the 5*0 (3) invariant metric. If we introduce a densitised 
version of S^ by 

j^ap = '^e"^^'S%s9AB, (2.2) 

this leads to the Lagrangian density 

^ = R^p^f- (2.3) 



2+2 Ashtekar variables II 3 

By making a 2 + 2 decomposition of the connection according to 

r^ = r^^dx^ = A^dx' + B^Jx"". (2.4) 

and using the S'0(3) covariant derivative D this may be written as 
C = S/Afo + ti'B^, + B^D^ti' + B^D.tX' + R\tl^ + R^t^^ (2.5) 

In paper I we considered the phase space to be given by the twenty three variables 
Sj*^, S^*, ii\, s\ together with the connection terms Af and B^. These twenty 
three variables were required to satisfy thirteen constraints which left ten degrees of 
freedom: two spin and boost freedoms (self-dual and anti self-dual) and eight degrees 
of freedom for the double null metric (the 10 degrees of freedom for a general metric 
are reduced to 8 because of the two slicing conditions). However, as explained in the 
introduction, in this paper we will work directly with the densitised bivectors and use 
all the components of f^^ rather than use a mixture of the sigma variables and frame 
components as we did in paper II. There are a total of 18 independent components and 
these have to satisfy nine constraints. This means that we have only nine degrees of 
freedom; one less than in paper II. The reason for this is that the variables we are now 
using are all manifestly self-dual so that we need only consider the effect of the self-dual 
spin and boost transformations on the sigma variables. This is different from paper II 
where the frame variables ^\ were not invariant under anti self-dual transformations 
and this therefore had to be considered as a gauge freedom. However in both cases the 
full Hamiltonian analysis which involves splitting the constraints into first and second 
class constraints gives the same number of dynamical degrees of freedom. 

We now give the nine constraints. The first two constraints are found by expressing 
the results S'^ A S"^ = and S*-*- A S*^ = in terms of S^ . This gives 

eaf3'y5^i"^^2'^ = 0> (2-6a) 

e„^,5Si"^I]3^^ = 0. {2M) 

The next set of constraints are obtained by considering the expressions for S J* in terms 
of the frame variables. This gives the following four constraints: 

ei,f:^'>%'' = 0, (2.7a) 

ei,t^'%'^ = 0, {2.7b) 

S^oi = 0, (2.7c) 

Sg"^ = 0. {2.7d) 

Two further constraints are found by expressing /x''"o/i°o and fJ^^if^^i in terms of the sigma 
variables. The double null slicing conditions (see paper II equation (2.23)) then give 

eijf:^'%'' = 0, (2.8a) 

6,,E2°^E3°^' = 0. {2.8b) 

We may now obtain a simplified version of these constraints by noting that 
()2.fia^ . ()2.7a|l . ()2.75j) . ()2.8a| and ()2.8&j) can be combined with the requirement that 
the volume form be non-vanishing 

^2 ^3 ^ab^ij 7^ U, 



Ci 




C2 




Cs 


= s,°i = 0, 


c. 


- S3°^ = 0, 


c. 


= S,°2 = 0, 


c. 


= S,°3 = 0, 


C-j 


- (Sa^^)^ = 0, 


Cs 


- i^s'r = 0. 
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where they reduce to the requirement that either 

t^o^ = S^o^ = S3I2 = S3I3 = 0, (2.9) 

or 

1:302 = t^03 ^ f^n ^ f^u ^ Q (2.10) 

These two conditions are interchanged when one swaps the x" and x^ coordinates so 
there is no loss of generahty in choosing the former condition which coincides with the 
choice of shcing condition made in paper II. We therefore have the following constraints. 

0, (2.11) 

0, (2.12) 

(2.13) 

(2.14) 

(2.15) 

(2.16) 

(2.17) 

(2.18) 

Note that we have chosen to square the final two constraints as these correspond to 
the existence of cychc variables and as observed by Goldberg et al ||5j squaring such 
constraints considerably simplifies the canonical analysis. The final constraint comes 
from expressing the constraint for the conformal factor C (see paper II equation (3.10)) 
in terms of the sigma variables. This results in the constraint 

C<j = tab^X ^1 + ^ab^ij^2^^^ "^1 ^1 • (2.19) 

We may now consider our Lagrangian to be given by ()2.5p . where the variables are 
the components of E^ and the connection F"^, and use Lagrange multipliers to introduce 
the primary constraints given above. This results in a Lagrangian density given by 

C = tX'At, + t^'B^^, + < (DiS^ + D,t^) + R\tX' + R\,tl^ - X-C^, (2.20) 

where a sums from 1, ..., 9. It may be verified that making a variation of this Lagrangian 
with respect to the components of S^ leads to the Einstein equations while variation 
with respect to the connection leads to the structure equations. We will not show this 
explicitly here as the calculations are very similar to the derivation of these equations 
in the Hamiltonian description which we give below. 

3. Hamiltonian description 

The Lagrangian density is of the form C = p^q\ — Ti, and therefore we can see directly 
from ()2.20|) that the canonical variables are A^ and B^, and have the respective 
momenta Sj** and Sj*^. We can also simply read off the Hamiltonian density which is 
given by 
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where we have introduced the momenta P"^^,P"^^,Pa for the cychc variables 
E^^, S^*, B^. This results in additional primary constraints which have been 
introduced into the Hamiltonian using the Lagrangian multipliers ^a^,^a\^a- Note, 
that unlike paper II we will not use the 'shortcut' method of treating the cyclic variables 
as if they were multipliers and eliminate them from the canonical analysis, but will 
include all the variables in the constraint analysis. This is because in the present case it 
is not at all obvious that the corresponding constraints are automatically propagated. 
The canonical Poisson brackets are given by 

{Afix),t^\y)] =6^SiSix,y), (3.2a) 

{Bfix),t^\y)] =S^Six,y), (3.26) 

{B^ix),PB{y)} =5i5ix,y), (3.2c) 

{t2Hi),P^23iy)} =S^Six,y), (3.2d) 

[t2%x),P%{y)} =6^6}6{x,y). (3.2e) 

We have a total of twenty one primary constraints introduced into the Hamiltonian. 
Following the Dirac-Bergmann algorithm we propagate the primary constraints. The 
equations for Ci, C*2 and Cg are rather complicated so we do not give them explicitly 
however they define the multipliers ^2^^) ^3^^ ^^^ ^1^^ respectively. The remaining 
constraints give 

Cs =[D,t2' + 2r]\sB^f:^')5^, (3.3a) 

C, =(D,t2' + 2r^\^B^tS')6f, {3.3b) 



^ Oi 

(75,6 = S2 = -2r]%^B^tS' - t^'], + 2r]^B2BF^A + D,t^'\ (3.3c) 

^7,8 = S3'' = ^3'^ (3.3c?) 

From which we see that the equations C7 and Cs define the multipliers ^3^*. We also 
need to propagate the momenta conjugate to the cyclic variables. This gives 

Pa =D^tX' + Dit^\ (3.4a) 

Pis = Ris - Sf{^2'>^' + Ss^'A^ - t^'^X^) - b^t^^X^ - btt^^\\ {3Ab) 
P\^ = R\^ - 6^{t^'^X' + SgO^A^ - Si°3a9) - 6^{t^'^X' - SgO^A^) 

- (5^(I]i03a2 _ s^o^A^ + 2S3i2A^), (3.4c) 
p^3 = i?^^3 + ^^^(S^o^Ai + SgO^A^ - Si°2a9) + 6^{t-,'^X' - Sg^^A^) 

+ 5^(I]i°2a2 - S2°2a9 + 21:31=' A^). (3.4d) 

Equations ()3.46|) - ()3.4dj) define the multipliers A^, A^ and A^ through the equations 
AiSi°i = R^s, A^S^oi ^ ^3^^ ^^^ ^9f;^oi ^ -R\^. This leaves thirteen secondary 
constraints given by 

Cs = f ASi* + 2r]''^sB^tS') 6^ = 0, (3.5a) 



varial 


Aes II 


6 


C, = 


(ASi' + 2r/^^s5o^S^"i) ^3^ = 0, 


{3.5b) 


X. Oi 
S2 = 


-2v'lAB^f^S' - S2':i + 2v%Bf^2' + ^.S^'^^' = 0, 


(3.5c) 


P\. - 


R\A'' - R\A'' - R%^2' - R%^s' ^ 0, 


{3.5d) 


n- 


R\A'' - R%^i' - R\,^s' = 0, 


(3.5e) 


p3 ^ 

-f^ li ~ 


R'u^r - R%^i' - R\j^2' = 0, 


(3.5i) 


Pa = 


Dit^^ + D,t^' = 0. 


(3.5^) 



The Dirac-Berginann algorithm then requires us to propagate these secondary 
constraints to see if any further tertiary constraints arise. Before this is done we rewrite 
the secondary constraint ( 3.5j| ) by looking at its components. This gives 



^^''P\ = ^s'\R\^i'' - R%^i' - R\^2'), (3.6a) 

Si^^p^H = f:^'\R\±^'' - R%f:,'' - R\f^2'')- (3.66) 

As we show in section 4 the five equations (j3.5dj) . ()3.5e|l and ()3.6a|l define five of the 
Einstein equations and are therefore conserved by the Bianchi identity. Propagation 
of p.5gD is identically zero for A = 0, and defines the multipliers A^ and A^ for the 
remaining values of A. Equations ()3.5aj) and ()3.56|) define the multipliers ^^ and ^^ when 
propagated. Propagation of ()3.5cjl define the multipliers ^2^\ ^^^d finally propagation of 
(j3.6fejl leads to a multiplier equation giving the multipliers .^j^^*. Therefore no additional 
constraints arise through the propagation of the secondary constraints. 

Having dealt with the constraints we next calculate the equations of motion. 
Evolving A^ gives 

Af = D,B^ - 6^e,,{t,'^X' + t,'^ X' - t,'^ X^) 

-8^e,,{tJ^X'~tJ^X^) (3.7) 

— 03 ejj(S-^ A — S2 A ), 

which results in 

A] = D.Bl - [R%t,'^ + R%t,'^ + R\^t,'^) /t,'\ 

if = D^Bl - {R%t,'^ + R\t,'^)/t,'\ 
if = a^Bl - {R%f:,'^ + R\t^'=)/t-,'\ 
The remaining equations of motion are given by 

^ Oi _ _ _ .. 

Sa = 2r^^BA<Sc°^ - DiSi^ - D,J:i\ 
B^ = D^B^ + dtit^^^'X^ + ti^X"" - Si^^A^) 
+5^{t-i'X^ + A=^) + dtiti'X^ + A^). 



(3. 


8 a) 


(3. 


86) 


(3. 


8 c) 


(3. 


Sd) 


(3. 


8e) 


(3 


8i) 


(3. 


H 
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4. Einstein Equations 

Since we have now obtained all the constraints and equations of motion we are now 
in a position to derive the Einstein equations. We use five of the secondary equations 
(j3.5d| I3.5el and l3.6ajl . together with the expression for the Einstein tensor in terms of the 
self-dual curvature ( see equation (5.1) of paper II) to obtain the following five Einstein 
equations. 

Oi' 



p2 y ^J - f 


p3 y Oj _ f 


pi y ^J - f 

^ lj^2 - (^ 


pi y oj r^ ( 


p2 y Oj ^ [ 



Oi 



732 v' 01 I d2 v Oi I d1 v 

K yi^i + K ijh-^ + K -jh^ 

p3 v-i 01 I p3 v-i Oi I pi v 
pi v 01 I pi v-1 Oi 1^ p3 v-i ( 



pi v 01 I pi V' Oi I p3 v' Oi 

K ^jh-^ + K .^2.1 + K .^.2.3 



p2 v-i 01 I p2 y-i Oi I pi V' 



Oi 



t,'^ ^ 
t,'^ ^ 
t,'^ ^ 



V Oi _ 

^3 ~ 







S °^' ^ 



G°o^0,(4.1a) 
G\^0,iA.lb) 
G\ ^ 0,(4.1c) 
G^a ^ 0,(4. 1(/) 
G^3^0.(4.1e) 



Ij^l "T -"- ij^l ">" ^^ ij^3 

The remaining Einstein equations are obtained from the equations of motion from which 
we obtain 



A^s^" ^ - {r\A'' + R\tr'' + R^.t^'') t^'' ^ ^^ G^ ^ 0, 
A]t^'^ ^ - {R\t-,'' + R\t^'^ + R^t^'^) t^'' ^0^^G\^0, 
if s,^^ ^ - {r%A'' + R\^2' + R%^,'') ^2' ^0^G\^0, 
Ajt,'^ + Bit,'' ^ - (^R\t,'' + R\t,'^ + R%t,'^) E,°^ 

pi V 01 pi V 23 p2 fi 23 p3 fi 23A v 01 ^ n > ^ n^ ~ 

n 01^1 -fl 23^1 ^ 23^2 ^ 23^3 I ^1 ~ ^ "^ ^ "-^ 3 ~ 



0. 



(4.2a) 
(4.26) 
(4.2c) 
(4.2(i) 

(4.2 e) 



We have therefore shown that we can derive all the Einstein equations from the 
constraint and evolution equations. The structure equations are obtained in exactly 
the same way as in paper II so we do not repeat the calculation here. The final stage is 
to ascertain which constraints are first class and then calculate the algebra of first class 
constraints. 



5. First class constraints 

We now move on to the next stage of the Dirac-Bergmann algorithm and calculate the 
first class constraints. This is harder than in paper II because we do not have available 
such an obvious geometric interpretation of the constraints. Fortunately some of the 
secondary constraints are the same as the corresponding secondary equations obtained 
in paper I and we can therefore adapt them by adding appropriate combinations of 
the other constraints in the same way as we did in paper II to obtain four first class 
constraints. However because we are not using the 'shortcut' method but including the 
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cyclic variables in our analysis we have to add extra terms to the three constraints ipi, 
and ipp to ensure they are first class. This results in the following constraints. 



+B^,Pa - tl^p^ 



23,1 



^aP\i + (Si^Pt).i ^ 0, 



i^V--=Bt^A+^t^A 



{Aft2' 



__ I A Ay 01^ 



+b^,Pa - tl'p% - tx^p\^ + {tx^p%] 



(5.1a) 



(5.16) 



We also need to modify the Gauss constraint to take account of the extra cyclic variables. 
This results in 



^1 := DiS^oi + ASi°^ + 2r,%^P^,^t^^ + 2ri''^^P%tS' - 2v%B^Pc ^ 0. (5.2) 



Because we are including the cyclic variables in the analysis it turns out that we also 
obtain an additional two first class constraints which are given by 



Pi 



23 



Pi 



0, 
0. 



(5.3a) 
(5.36) 



We can show that there are no further first class constraints so that we have a total of 
six first class given by ()5.16|) - ()5.36|) . The remaining constraints are second class which 
for completeness we list below; 



pA 

-* li 



P- 
P- 



a 



1,...,9 



23 



0, 
0, 
0, 



23 ~ "' 
^2 = 0, 

^3 = 0, 



01 
'A 



DiJ: 



01 
A 



D,t2' 



D,t2' 



0, 
0, 



D,tX' + 2r,''j,^B^tS')5^ 



D,t}'^ 



^ri'^ABB^tS'^St 



0, 
0, 



R'\,^r + R%T.^'']5, 



TjA V^ 01 , TjA ST' I 



^7 2A^0 ^C 






0. 



(5.4a) 
(5.46) 
(5.4c) 
(5.4d) 
(5.4e) 
(5.4/) 
(5.4^) 

(5.4/i) 

(5.4^) 

(5.4j) 

(5.4A;) 

(5.40 

(5.4m) 
(5.4n) 



Because we are not using the 'shortcut' method but working with a phase space which 
includes all the cyclic variables we have a rather large phase space consisting of 42 
functions. These functions are subject to 6 first class and 28 second class constraints 
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leaving 2 dynamical degrees of freedom per hypersurface point as is appropriate on a 
null surface 0, [7j, [HI, P. 

Having obtained all the first class constraints we calculate the first class algebra by 
computing the Poisson brackets of the constraints with each other. To do this we first 
obtain smooth versions of the constraints by smearing them with smooth functions. We 
start by considering ipi. Let F be a vector field with components tangent to {5*} and 
define a smeared version of the tpi constraint by 



^(F) 



FVi d^x. 



(5.5) 



We next consider ipi. This time we let F be a vector field with components tangent to 
the null generators of Sq and defined a smeared version of the ipi constraint by 



^(F) = / F^iPid^x. (5.6) 

In a similar way we define a smeared versions of the remaining constraints by 

G{f)= [fGid'x, (5.7) 



Pif) 

n(/) 



fPid'x, 
fP\, d''^ 



(5.8) 
(5.9) 



where / is a scalar field and / a scalar density. We now obtain the first class algebra by 
calculating the Poisson brackets of all the first class constraints with each other. Below 
we show only those terms that are not strongly zero 

(5.10a) 
(5.106) 
(5.10c) 
(5.10c?) 
(5.10e) 
(5.10/) 
(5.10^) 
(5.10/i) 

(5.10i) 

As in paper II we have chosen to keep Vi ^"^d ipi separate to illustrate the 2+2 
structure of the constraint algebra. However they may be combined to give ipA, where 
ii'A) = {i'1,4'2,4'3)- This may be smeared with a general vector field F on Sq to give 

^(F)= [ F^ijAd^x. (5.11) 



^{R),^{si^ 


> = ^{Cj,S) 


^{R),^{S)\ 


= ^(^rS), 


^iR),^iS)} 


= ^(^rS), 


^iR),Gis)} 


= G{Cf,s), 


^iR),Gis)} 


= G{C^s), 


§(/?), H(.~)} 


= n(/:^s), 


^(i?),H(5)} 


= n(£^s), 


^(i?),P(S)} 


= Pi^R-s), 


^{R),P{~s)} 


= Pi^RS). 
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The constraint algebra then has the more compact form 

{^iR),^iS)} = ^iCRS), (5.12a) 

{^{R),G{s)} = G{Cns), (5.126) 

{vi/(i?),n(5)} = n(/:R5), (5.12c) 

{■^iR),Pis)} = PiCRs). (5.12d) 

Having calculated the first class algebra, we give a geometrical interpretation of the 
first class constraints by calculating the infinitesimal transformations the constraints 
generate. However we will not consider the constraints Pi = and P\^ = (which 
come from the cyclic variables) as they just indicate the gauge freedom to choose the 
variables Bq and S^^^. We start with the constraint ipi. Calculating the Poisson brackets 
of this with the other variables gives 

Mi = {Af,^iG)] = CgA^^, (5.13a) 

5B^ = [B^,^iG)] = CcBi, (5.136) 

5B^ = {i?o^,^(G)} = CcBi, (5.13c) 

6t2' = {si^^(G)} = Cg^2^ {5.13d) 

5t2' = |si\^(G)) = £GSi\ (5.13e) 



From this and we see that, as in paper II, ipi generates diffeomorphisms in the two 
surface {S}. 

We next consider the constraint ipi. 

5Ai =\Af,^iG)}=CGAi, (5.14a) 



5B^ = [B^,^iG)] = CoBi, (5.146) 

5B^ = {5o^,^(G)} = CgB^, (5.14c) 

6t2' = {si^^(G)} = Cg^2^ (5.14d) 

5t2' = (si\^(G)| = Cat^'. (5.14e) 



From this we see that this constraint generates the diffeomorphism along the null 
generators on Sq. 

Finally we look at the commutators with the modified Gauss constraint Qi. These 
give 



6Af 


= {Af,G{g)} = 


-g,6t - 2gA^6^ + 2gAf6^, 


(5.15a) 


6B^ 


= {Bf,G{9)} = 


-g,^6^-2gB',6^ + 2gB',6^, 


(5.156) 


6B^ 


= {B^,Gig)} = 


-2gB^X + ^gB^^t 


(5.15c) 


6t2' 


= {^r,Gi9)} = 


- -2gtiH^ + 2gtiHt 


ib.lbd) 


S^A 


= h^\G{g)} = 


-2gt^'5^ + 2gt^'5-^. 


(5.15e) 
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Comparing these infinitesimal transformations with the effect of the self-dual spin and 
boost transformations (see equation (7.15) of paper II) we see that Qi again generates 
the self-dual spin and boost transformations for these variables. 

6. Conclusion 

In this paper we have obtained a 2+2 double null Hamiltonian description of General 
Relativity using only the complex 5*0(3) connection and the components of the complex 
densitised self-dual bivectors S^. We have carried out the general canonical analysis of 
this system and obtained the first class constraint algebra entirely in terms of the self- 
dual variables. All the first class constraints have a simple geometrical interpretation. 
The constraint ipi generates the diffeomorphisms in the two surface {S}, and the 
constraint ipi generates the diffeomorphisms along the null generators of Sq, so that 
these constraints can be combined to give ipA which is just the momentum constraint 
in this description. The modified Gauss constraint Qi generates the self-dual spin and 
boost transformations, while the final two constraints correspond to cyclic variables and 
express the freedom to freely choose the corresponding canonical variables. The choice 
of variables we have made makes it possible to compare the results of using a double 
null evolution with both the standard 3+1 approach using Ashtekar variables ITU] and 
the null approach of Goldberg et a/ [S]. In particular one can see how the form of ipA is 
very similar to that of the adapted first class version of the momentum constraint 



^ijK>^''-AlW,(U^:.)=0 (6.1) 

which is used in the standard Ashtekar approach. 

The next step of the canonical quantisation process is to eliminate the second class 
constraints by replacing the Poisson brackets by Dirac brackets [TTj. These are modified 
versions of the Poisson brackets such that the Dirac bracket between any of the second 
class constraints and any other variable vanishes identically. The details of the analysis of 
the second class constraints and the construction of the Dirac brackets will be presented 
elsewhere, but we outline the general procedure below. We start by introducing a 
vector K whose components Kj, I = 1, ... 28 are the second class constraints given by 
equation ()5.4aj) - ()5.4?ij) . We then calculate the Poisson bracket matrix C for the second 
class constraints whose components are given by 

Cjj = {Cj,Cj}. (6.2) 

Although in principle C is a large matrix, the form of the second class constraints mean 
that, as in paper II, many of the entries are identically zero. 
Having calculated C the Dirac brackets are then given by 

{F, G}d = {F, G} - 5^{F, Cj}C-''^{Ck, G}, (6.3) 

J,K 

where C~^ is the inverse of the matrix given by ()6.2|1 . Note that the invertibility of 
C is guaranteed by the independence of the second class constraints. Rather than 
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calculating the Dirac brackets an alternative procedure would be to impose a suitable 
gauge condition and use this to explicitly solve for the second class constraints and hence 
obtain a Hamiltonian which only depended on the true dynamical degrees of freedom 
(see e.g. Goroff and Schwartz [T2]) 

As well as working with entirely self-dual variables another important feature of the 
approach taken in this paper is the use of a double null evolution. This has the advantage 
over choosing a null hypersurface that the projection operation is well defined. It also 
has the advantage over the standard 3+1 approach that the Hamiltonian constraint 
is second class rather than first class. This has the important effect that the first 
class constraint algebra also forms a Lie algebra. As described above the next step of 
the standard canonical quantisation process would be to calculate the Dirac brackets. 
However an alternative approach would be to look at the corresponding loop variables 
on the null hypersurface. This might provide a link between the Ashtekar approach 
to loop quantum gravity (see for example ^H]) and the work of Iyer, Kozameh and 
Newman on formulations of Einstein's equations based on the holonomy of null surfaces 



d'Inverno, R. & Vickers, J. (1995)- '2+2 decomposition of Ashtekar variables'. Class. Quantum 

Grav. 12(3), 753-769. 
d'Inverno, R. A. , Lambert P. & Vickers, J. A. (2005). 'Hamiltonian analysis of the double null 
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